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Current fluctuations in non-equilibrium diffusive systems: an additivity principle 
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We formulate a simple additivity principle allowing to calculate the whole distribution of current 
fluctuations through a large one dimensional system in contact with two reservoirs at unequal 
densities from the knowledge of its first two cumulants. This distribution (which in general is non- 
Gaussian) satisfies the Gallavotti- Cohen symmetry and generalizes the one predicted recently for 
the symmetric simple exclusion process. The additivity principle can be used to study more complex 
diffusive networks including loops. 
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Understanding the fluctuations of the steady state cur- 
rent through a system in contact with two (or more) heat 
or particle reservoirs is one of the simplest and most fun- 
damental problems of non-equilibrium physics Q, Q, |^ . 
For quantum particles such as (weakly interacting) elec- 
trons which satisfy the Pauli principle, the whole dis- 
tribution (the full counting statistics) of the number of 
particles transferred between the two reservoirs during 
a long time interval is known and it can be calcu- 
lated by a number of theoretical approaches 1^ , 
ranging from the theory of random matrices [^Jg] to the 
Boltzmann-Langevin semiclassical description (9^. 

For systems of purely classical interacting particles 
P, 01 in contact with two reservoirs the theory is, to 
our knowledge, less developed. However, for a num- 
ber of stochastic models of classical interacting particles 
[lol im , the cumulants of the current fluctuations 

were found to coincide with those previously known of 
non-interacting quantum particles. It is of course an im- 
portant issue to know under what condition a classical 
particle system could present the same distribution of 
current as in the quantum case. 

For most theoretical approaches developped in the 
quantum or in the classical description, the calculation of 
the cumulants becomes harder and harder as the degree 
of the cumulants increases. The goal of the present let- 
ter is to show that for classical stochastic models, if one 
postulates a simple additivity principle for the current 
fluctuations, the whole distribution of current fluctua- 
tions can be calculated from the knowledge of the first 
two cumulants of the current. 

We consider here a one dimensional diffusive open sys- 
tem of length N (with TV large) in contact, at its two 
ends, with two reservoirs of particles at densities pa and 
Pb. In the bulk, the system evolves under some conserva- 
tive stochastic dynamics and, at the boundaries, particles 
are created or annihilated to match the densities of the 



reservoirs. 

Let Qt be the integrated current up to time t, i.e. the 
number of particles which went through the system dur- 
ing time t. For large N, we shall see that the whole 
distribution of the fluctuations of Qt depends only on 
two macroscopic parameters D(p) and cr{p) deflned as 
follows: Suppose that for Pa — P + and pb — p with 
Ap small, we know that in the steady state Fick's law 
holds 



{Qt 



(1) 



Suppose that for Pa = Pb = P (in which case (Qt) — 0), 
we also know that for large t the variance is 



{Q' 



(2) 



The main result of the present paper is that, using a 
simple additivity principle H1UI11|) , we can predict all the 
cumulants of Qt for arbitrary pa and pi,. If we define the 
integrals /„ by 



D{p) a(p)"-i dp, 



the first cumulants of Qt are given by 



(Qt) 1 , (g?) - {Qtf 



N h 



1 3(13/1 - /|) 



N 



t 



1 3(5/4/f - 14/1/2/3 + 9/|) 



N 



(3) 



(4) 



(5) 



The case pa = pt can be obtained by letting pa tend to 
Pb in the previous expressions. 
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More generally, all the higher cumulants can be ob- 
tained from the knowledge of /ijv which characterizes the 
large t growth of the generating function of Qt 



fiN{X,Pa,Pb) = lim t-Hii{e^Q') 

t—>-oo 



(6) 



We are going to show that, for large N, takes the 
following parametric form 



PN{\Pa,Pb) = --^ 



D{p) dp 



H 2 



/p, ^l + 2Ka{p) 
where K = if(A, pa, pb) is the solution of 



+ 2Ka{p) 



(7) 



(8) 



As pN - iHQt) + >^''{Qt)c/2 + X^Q!)c/Q + one 
simply needs to expand {T)) and © in powers of K and 
to eliminate K to obtain fi^ as a power series of A and 
the cumulants such as (jSEJ- 

Note that and ((HJ are only valid for pa ^ Pb and 
in the range of values of A where K is large enough for 
the argument of the square root in the integrants not to 
vanish. We checked that they can also be analytically 
continued to cover the other ranges of A and the case 
Pa = Pb- 

Our derivation of and ^ is based on an additivity 
principle that we are going to formulate now. The prob- 
ability Pn{<1-, Pa, Pb, t) of observing an integrated current 
Qt = qt is exponential in t for large t 

PNiq,Pa,Pb,t) exp[t Fn {q, Pa, Pb, conta,cts)] , (9) 

where Fjv((j, Pa, Pb, contacts) depends on the length N of 
the system, on q, on the densities pa and pb in the two 
reservoirs, and on the nature of the contacts of the sys- 
tem with the two reservoirs. {Fn is negative and vanishes 
only when q takes its most likely value {Qt)/t). When 
N is large and q is of order 1/N, the effect of the con- 
tacts becomes negligible and asymptotically FN{q, pa, Pb) 
depends only on q,pa,Pb, on the length N and on the 
bulk properties of the system. We then assume that, for 
large N and q of order the large deviation function 
Fniq, Pa, Pb) satisfies the following additivity principle: 

FN+N'{q,Pa,Pb) ^ max{i^Ar((j, Pq,p) + FN'{q,p,Pb)} ■ 

(10) 

This property simply means that the two subsystems are 
independent, except that they try to ajust the density p 
at their contact to maximize the following product 

PN+N'{q, Pa, Pb,t) ma:>c[PN{q, Pa, P,t) Pn'iq, P, Pb,t)] ■ 
p 

We make also the following scaling hypothesis 

FN{q,Pa,Pb)^N~^ G{Nq,pa,Pb). (11) 



This hypothesis, which is valid in particular for the sym- 
metric simple exclusion process, means that pN defined 
by © is of order 1/N for large A'' (see [isjl. 

If we write N = {N -\- N')x, i.e. we split a system of 
macroscopic unit length into two parts of lengths x and 
1- X, then (|10lll|l lead to 



G{q,Pa,Pb) 



max 
p 



G{qx,pa,p) , G{q{l - x),p,pb) 



12) 



If wc keep dividing the system into smaller and smaller 
pieces and we use that for a piece of small (macroscopic) 
size Ax (i.e. of NAx sites) one has (P I2I10I11|I 



— G{qAx,p,p + Ap) 
Ax 



[qAx + D{p) ApY 



2a{p)Ax 
one finds a variational form for G 

"1 [q + D{p{x))p'{x)y 



G{q,pa,Pb) = -min 

p(x) 



2a{p{x)) 



dx 



(13) 



(14) 



where the minimum is over all the functions p(x) with 
boundary conditions p{0) — pa and p{l) = pb- 
The optimal p{x) in H14|l satisfies 

l'a'{p)-c'{p) (^gy_2c(p)g=0, 

where a{p) = {2a{p))~^ and c{p) = D'^{p) a{p). Multi- 
plying the above equation by dp(x)/dx, one obtains after 
one integration 

D'ip)(^^J = q\l + 2Ka{p)), (15) 

where i^T is a constant of integration. 

To proceed further one needs to determine the sign of 
The simplest case is when p{x) is monotone, and this 
happens when q is close enough to its average value for 
Pa 7^ Pb (this corresponds to values of K small enough 
for the right hand side of 1)15(1 not to vanish). The in- 
vestigation of this regime is enough to determine all the 
cumulants. If for example pa > Pb, the optimal p{x) is 
decreasing for small K 



dp 

dx 



Dip) 



Vl + 2Ka{p) , 



and this leads to the following expression for G: 



G 



"° Djp) 
. ^ip) 



1 + Ka{p) 
Vl + 2Ka{p) 



dp. 



where the constant K is determined by: 

q=rdp^£^. 

Jp, ^1 + 2K<j{p) 



(16) 



(17) 



(18) 
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One can then show that 

dG G ^ D{p) 

oq q Jp, cr{p) 



1 - 



VI + 2Ka{p) 



where the derivative is taken keeping pa and pb fixed, and 
using the fact that /ijy = N'^ maXq[Xq+G{q, Pa, Pb)], one 
obtains l(715|l . 

When the optimal p{x) is no longer monotone, i.e. K 
is negative enough for the right hand side of (|15|l to van- 
ish, the expressions (|7I8I17I18() of pN,^,G,q are modi- 
fied. We checked that their new expressions are simply 
the analytic continuations of (|7I8I17I18(I . 

In general when the system is in equilibrium (pa ~ 
Pb ~ p) the fluctuations given by (|14|l are non Gaussian. 
However when Pa = Pb — P* where p* is the density 
for which cr(p) is maximum, the optimal p{x) in \\A\ 
satisfies p'{x) = and the fluctuations become Gaussian 
{G{q,p*,p*) = -q'^ /(2a(p*))) in agreement with the con- 
jecture made in for a specific model, the symmetric 
simple exclusion process. 

It is also easy to check from (|14() that the optimal pro- 
file p{x) is the same for q and —q. This implies that 

G{-q, Pa, Pb) ^ G{q, Pa, Pb) - 2q [ ^^r-^dp 



Jp, <P) 

which is the Gallavotti-Cohen relation 0, Q, 0| 



a linear system with a length given by Kirchoff's law for 
the addition of resistors. 

We consider now two specific examples of stochastic 
dynamics on a Id lattice, the symmetric simple exclusion 
process (SSEP) and the zero range process (ZRP). The 
number of particles at site i G {0, A''} is denoted by rji. 

In the SSEP, each site is is either empty or occupied by 
a single particle (?7i = or 1) and each particle attempts 
to jump to its right or to its left at rate 1 if there is no 
other particle at the corresponding neighboring site. Let 
Q\ be the integrated current through bond {i, i-\-l) during 
time t. As in the steady state the integrated current is 
independent of the bond and dt{Q\) — {rji — jyi+i), 

NdtiQt) =^{'ni- Vi+i) = Pa - Pb, 

i 

thus D{p) = 1. As N^dtiiQt)^) = we 
write 



+ - V^+l)) + {V^+l{'^ - V^)) 

The first term simplifies 





£l 


Pa ' 






U 




Pb 





FIG. 1; The system connecting the two reservoirs contains a 
loop with two arms of unequal lengths. 



Consider now a system composed of 4 parts as in Fig- 
ure n The left reservoir is connected to C by a chain 
of length Nxi. Between C and D there is a loop made 
of two chains in parallel of lengths Nyi and Ny2 and D 
is connected to the right reservoir by a chain of length 
Nx2- According to the additivity principle, one should 
have 



Gioop(g, Pa, Pb) = maxp^^ 



Pd,q' 



+ 



G{q'vi,Pa,Pd) I G{(q-q')v2;Pa,Pd) , G{qX2,Pd,Pb) 
VI y2 X2 



The optimum is achieved when q'{yi +2/2)= <Z2/2, thus 

Gloop(g, Pa, Pb) = G{qu, Pa, Pb)/u , 

with u = xi + (j/f ^ -I- y^^)~"^ + X2- So the current fluc- 
tuations for the system with a loop are the same as for 



and it vanishes for p = pa = Pb- For pa = Pb, the sta- 
tionnary measure is product so that a{p) = 2p(l — p) 
according to lj2Jl. The cumulants derived in [l^l as well 
as the expression conjectured for p^ 

AtJv(A) = -iV"i[sin"i(V^)]2 , fortj<0. (19) 

where uj — {\ — e^^){e^pa — Pb ~ {^^ ~ ^)PaPb) can be 
recovered from (|7I8|I . This can be seen by noticing that 
the optimal profile solution of (|16|l is 

_ 1 / sin {2{ea + {9b -9a)x)) \ 
P^''>-2['+ sin(2/) j' 

where the parameters /, 9a , 9b are fixed hy K = 
tan^(2/), p(0) = pa,p{l) = Pb- In terms of these pa- 
rameters, A and p^ take the form 



A = loE 



cos(/ + 9a) cos(/ - 9b) 



cosif- 9a) cos{f + 9b) 



PN = -i9a~9b)'' 



For the ZRP the number of particles on each site can 
be arbitrary and the jump rate ^(jji) from site i to each 
of its neighbors is an increasing function of the number 
of particles rji at this site. We choose $(0) =0. We have 

J2dt{Ql) = J2{^v^) - Hm+i)) = *(pa) - ^(Pb) , 
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where the expectation of $ under the stationary measure 
at density p is denoted by ^'(p). We also have 

N^mQtf) = Y.^Ql{Hv^) - <f fei))> + 2^(<i>(r;,)) • 

hi i 

As for the SSEP the 1st term in the rhs of the above 
equation vanishes when pa = Pb- Wc finally obtain 
D{p) = a'{p)/2 and <7{p) = 21'(p) according to ifTI^ . 
Therefore 

Mjv(A) - (1 - e-^){e^a{pa) - a{pb))/2N . 

This generalizes the case of non-interacting particles for 
which <7{p) = 2p. The optimal profile is obtained by the 
change of variables 



a{p{x)) 



{0b - ea)xf - 1 

2K 



where 9a, Ob are fixed by p(0) = pa and — pb- In 
particular, the expression of pn follows from 



A = log 



1 



1 



Mw(A) 



9a - Obf 



The additivity principle H9I10|I formulated here and 
its variational expression (|14|l can be derived (work in 
progress) f rom the hydrodynamic large deviation theory 
llllliallJ- This theory was extended recently by Bertini 
et al. [13 who could calculate the density large deviation 
functional of the steady state as the optimal cost per unit 
time for a space/time density fluctuation. For diffusive 
systems, the exponential cost of observing an atypical 
space / time density profile over a time t can be estimated 
by a functional depending only on D(p),a(p) and on the 
density {pix, s)}^^[o,i],o<s<t (see eg [13 El El)- The 
optimal strategy to create a fluctuation of the current 
Qt = qt over a very long time t is to create a fixed den- 
sity profile p{x) in order to facilitate the deviation of the 
current and (|14|l can be understood as the cost for main- 
taining this atypical density profile. The optimal profile 
controling here the current fluctuations is time indepen- 
dent, in contrast to the one which controls the steady 
state density fluctuations that Bertini et al. E3 had to 
calculate. This is why our task here was easier and the 
additivity principle (|10|l is simpler than the one obtained 
in for the steady state fluctuations of the density. 

It would be interesting to see whether the Bertini et 
al. macroscopic fluctuation theory satisfies a generalized 
additivity principle for time-dependent densities in the 
reservoirs. Other interesting extensions of the present 
work include the study of the effect of asymmetry in the 
bulk dynamics (i.e. of a field which favors jumps of parti- 
cles from left to right) [H [H [H H Ei or the analysis 
of more complex networks, in particular of systems in 
contact with three or more reservoirs [2^ 113, HO] ■ 

Of course, a very challenging issue would be to see 
whether the additivity principle could be valid for some 



mechanical systems satisfying (|1I2|) without the need of 
an intrinsic source of noise as the stochastic systems con- 
sidered here. 
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